Dipole-dipole instability of atom clouds in a far-detuned optical dipole trap 
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The effect of the dipole-dipofe interaction on the far-off-resonance optical dipole trapping scheme 
is calculated by a mean-field approach. The trapping laser field polarizes the atoms and the accom- 
panying dipole-dipole energy shift deepens the attractive potential minimum in a pancake-shaped 
cloud. At high density the thermal motion cannot stabilize the gas against self-contraction and 
an instability occurs. We calculate the boundary of the stable and unstable equilibrium regions 
on a two-dimensional phase diagram of the atom number and the ratio of the trap depth to the 
temperature. We discuss the limitations imposed by the dipole-dipole instability on the parameters 
needed to reach Bose-Einstein condensation in an optical dipole trap. 



S 1 
i 

a 

o : 
^— > 

o ■ 

CZ3 , 
>v 

43 ■ 

(N : 
> ■ 
o ■ 

00 1 
(N . 



o . 

o : 

43 : 

Or 



- 1—1 

X 



I. INTRODUCTION 

The optical dipole trap provides a conservative and 
tightly confining trapping potential for polarizable par- 
ticles. Owing to its simplicity it has become a standard 
tool for manipulating neutral atoms Today, various 
dipole trap based systems serve the study of many-body 
problems where the atom-atom interaction is an essen- 
tial ingredient. In this paper we consider far-off-resonant 
trapping of polarizable particles interacting radiatively 
through the electromagnetic vacuum field. 

In dense samples of cold atoms, two-atom processes 
strongly influence the trapping. For example, col- 
lisional processes (photo-association, hyperfine ground 
state changing collisions, etc.) are known to result in 
trap losses and they can limit the maximum achievable 
density. Collisions usually depend heavily on the internal 
electronic structure of the species. For interatomic dis- 
tances in the range of the optical wavelength, the atom- 
atom interaction is dominated by the radiative electro- 
magnetic coupling. In this case one can distinguish two 
limits: the photon scattering is dominated (i) by spon- 
taneous emission, e.g. in magneto-optical traps (MOT), 
and (ii) by stimulated emission, which occurs in far-off- 
resonant dipole traps (FORT). Since the MOT enabled 
us first to capture dense atom clouds from vapour, it was 
first analyzed with respect to many-body effects 0, 0|. 
In a MOT, the cooling laser is quasi-resonant with the 
atoms and the sample forms an optically thick medium. 
The depletion of the laser beams together with the mul- 
tiple spontaneous scattering of resonant photons within 
the sample can lead to instability [2J and extra heating 
d, At high densities, the reabsorption of photons in- 
fluences the laser cooling itself @, 0]. The combination 
of these effects can result in a highly nonlinear collective 
dynamics of the atom cloud in a MOT 0, Q . 

The FORT operates at an extremely low spontaneous 
scattering rate, thus the effects of multiple spontaneous 
photon scattering are strongly suppressed. The mecha- 
nism of trapping relies on the process of absorption and 
stimulated emission of laser photons. This process polar- 
izes the particles, hence the dipole force is accompanied 
by the dipole-dipole coupling between atoms. 



In previous works a phenomenological term propor- 
tional to the square of the atom density was introduced 
to describe the atom-atom interaction, e.g., the effect of 
collisions on the loading of a FORT [H], [HI ■ This domi- 
nates in strongly localized traps, collisional blockade can 
prevent us from confining even two atoms in an extremely 
tiny FORT [l2| • This paper is devoted to studying, based 
on a microscopic theory, the radiative interaction be- 
tween atoms in a FORT. We neglect the effects of the 
Coulomb and exchange interaction terms that would re- 
quire the use of molecular potential surfaces [l3]. We 
explore the limits of the maximum achievable density in 
a FORT, which arise solely from the dipole-dipole energy 
shift. 

Tuning the laser frequency wl very far below the 
atomic resonance lua makes it possible to eliminate the 
spontaneous photon scattering at a constant depth of 
the trap potential, since the scattering rate and the trap 
depth scale differently with the detuning Aa — lu-^ — da- 
The former is proportional to 27$7 2 / A\ while the latter 
to fi 2 /A J 4, where f2 is the Rabi frequency of the atom- 
laser coupling, 7 is the atomic linewidth (HWHM) p|. 
With temperatures well below the potential depth, i.e., 
H£1 2 /Aa S> ksT, the particles are captured in the trap 
for long times. The dipole-dipole energy shift scales simi- 
larly to the spontaneous scattering rate, however, it is en- 
hanced by the atom density. Thus it may become signifi- 
cant as compared to the dipole trap depth when the atom 
number TV in the trapped cloud obeys N^/Aa > 0.1. 

High atom densities can be an issue for many kinds of 
experiments with FORT's. One example is the attempt 
to achieve Bose-Einstein condensation with alkali gases in 
optical rather than magnetic traps [HI, EH • Another ex- 
ample is connected to the so-called cavity cooling method 
which is suitable to complement the deep conservative 
potential of FORT by a damping force induced by the 
field of a high-finesse optical resonator [l6|, [jjj • Various 
cavity cooling setups have been realized: (i) the FORT 
field can be along the cavity axis and detuned from the 
near-resonant cavity field by several free spectral ranges 
[la UM i or the FORT lasers can be perpendicular to the 
cavity axis with photons scattered into the cavity respon- 
sible for cooling [2(J, HH ■ The simultaneous cooling and 
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trapping gives rise to very long trapping times for single 
atoms. Cavity cooling can be applied to many atoms as 
well, and should work efficiently until they get localized 
at the antinodes of the red detuned trap field [22|, HH, H3 • 

The stability of an atom cloud against dipolc-dipole 
coupling is an issue also in the ultracold tempera- 
ture regime where the atoms form a degenerate Bose- 
condensate. The corresponding anisotropic potential can 
be included in the Gross-Pitaevki equation [2||. There 
is a stable solution for the atomic mean-field wavefunc- 
tion depending on the trap aspect ratio [26j], or on the 
scattering length [U [28| ■ The dipolc-dipole interaction 
gives rise to density modulations [29j |. solitons [30(. In 
these works the atoms possess a permanent dipole mo- 
ment. Recently, the effects of the magnetic dipole-dipole 
coupling has been observed in a chromium condensate 
[3l| . The laser induced dipole-dipole coupling at low 
temperatures has been treated in [321 ] and was shown 
to lead to nonlinear atom optical effects. In this paper, 
Bose-Einstein condensation is concerned inasmuch as we 
study whether the dipole-dipole interaction excludes the 
possibility of forming a condensate in an optical trap. 

The paper is organized as follows. We start by reca- 
pitulating the formulae of the dipole-dipole interaction of 
atoms at a microscopic level in Sec. HU Then, in Sec. Mil 
we take the far-off-resonance limit and keep only the 
leading order terms in the small parameter 7/A^. The 
basic equations of the mean-field model are introduced 
in Sec. IIV1 These can be expressed in terms of dimen- 
sionless parameters, the scaled intensity and the scaled 
atom number, which define the universality of this model. 
We present numerical results of the stability analysis in 
Sec. [V] There we show the parameter regimes of sta- 
bility and instability on phase diagrams. In Sec. I VII the 
relation of the dipole-dipole instability limit to the limit 
of quantum degeneracy (Bosc-Einstcin condensation) is 
discussed. We conclude in Sec. IVII1 We put in an Ap- 
pendix the derivation of the density-enhanced terms by 
a systematic expansion, leading to a Hamiltonian which 
accounts for the two-body interaction to all order. 



II. DIPOLE-DIPOLE INTERACTION 

The theory of radiative atom-atom interaction in the 
presence of a drivin g la ser field is described in several 
papers, e.g. Refs. 171 132. \33L [34j . we will mostly use the 
approach presented in [331 ] - 

We consider a number N of atoms interacting with a 
single Gaussian standing-wave laser field mode along the 
z direction which has a mode function 



f(r) = ecos(k z z) exp [-2ir 2 (x 2 + y 2 )/w 2 



(1) 



where e is the field polarization, and the mode is paraxial, 
k z 2n/w. The atomic transition frequency is coa, that 
of the laser field mode is lol , and the detuning is defined 
by — ojl — to a- The atom-mode interaction strength 
is described by the Rabi frequency Q in the position of 



maximum coupling. Wc assume an S <-> P transition 
with a degenerate manifold of excited states and keep the 
three-dimensional polarizability of the atoms. Actually 
the fixed field polarization selects two levels taking part 
in the dynamics. The atomic internal degree of freedom is 
described by the vectorial lowering operator a — J2 q ?qO' q 
with q = ±1 and q = corresponding to the circular and 
linear polarizations, respectively. The quantization axis 
will be defined in accordance with the choice of the field 
polarization e. 

The equation of motion for the density operator in the 
Markoff approximation reads 



p=^-[H,p]+£p. 
in 



(2) 



In a frame rotating at the laser frequency lol 1 the Hamil- 
tonian is 



N r 



2m 

N 



fry £ <4j9$ 

mi i. )<7n , (3) 



n,m=l 

where r„, p n and a n are the position, the momentum 
and the polarization of the nth atom (n = 1,...,N). 
Next to the single atom terms, i.e. kinetic energy, inter- 
nal energy, and atom-field coupling, the last term con- 
tains the induced dipole-dipole interaction energy of the 
atoms. Note that the natural linewidth 7 characterizes 
the strength of this interaction. The tensor (3 depends 
on the coordinate difference R m n = 4 — r n of the inter- 
acting pairs of atoms. 

The dipole-dipole interaction is mediated by the broad- 
band vacuum, and is therefore accompanied by incoher- 
ent evolution. This is represented by additional terms in 
the Liouville operator responsible for the dissipation: 

N 



c p = -i^Z p}- 

71=1 

[ d 2 u N q (u)ale- lkAS ^pe lkASP "a q n 

N 

"7 ({vtn a (Rrnn)Vn, P}~ 

n,m=l 



(4) 



where { , } denotes the anticommutator. The single 
atom terms include the spontaneous decay accompanied 
by momentum recoil. The tensor N(m) = (1 — u ou), 
and its diagonal elements N q (u) = e^N(u)e q are the an- 
gular momentum distribution of the spontaneous emis- 
sion from the g-state in the excited manifold. The dou- 
ble sum describes the loss effect due to the dipole-dipole 
coupling. 
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In free space the tensors a and (3 assume the following 
form : 



tion z. The atomic polarization is described by the op- 
erator er +1 , and the projection Eq. ([7]) gives 



\ 3 J /i £> £> ,sin kR m 

Ot{Kmn) — r \ (1 — n mn O K mn ) — — 

z I Kt\ ran 



(1 — 3R mn o R n 



COS kRrnn SUl kR n 



{kRn 



(kRrnri) 



(5) 



rtiii — ^ (1 3 COS Ornn) 



sin kR Tl 



cos kR n 



^cos 



(^mn) 2 (kRrnri) 



^ COS kRfnn (10) 
kRrnn 



, — ? 3 I , ± ^ . cos JxiI^jjifL 

P(Rmn) = — \ (1 — -Rmn ° fimn) TTL 

z Kt\, rari 



(1 — 3R mn o i? mn ) j 



fci?„ 



cos fci? r , 



\kRmri) ft-riin)' 



(6) 



The angle mn = /L(R mn ,z) is the angle between the 
distance vector of the two atoms and the axis z. The 
Fourier transform is 



(3(k) = lira 



(fc 2 + fc2) (jfc2 



•J-o 2^^ (k 2 - klf + 2if{k' 1 + k\) 



(11) 



where k = \k\ w fc z , i? m „ = |i?„ m |, and i? mn is a unit 
vector along the direction of R m n- 

The fixed field polarization selects the excited state 
and the atom reduces to a two-level system with a n = er n e 
(n = 1, • • • , -/V). The tensors a and (3 have to be pro- 
jected onto this particular polarization, 



a} n f}(R mn )a n = fliRmnWmCr-n 



(7) 



where (3{R mn ) = f3(R mn )e. We now evaluate this pro- 
jection in two cases: for linear polarization along x, and 
for circular one in the x-y plane. 



A. Linear polarization 

When the polarization of the beam is linear, e = x, the 
atomic quantization axis is taken in this direction. The 
atomic polarization is described by the operator ctq, and 
the projection given in Eq. © results m 



flmn = ~ (3 COS 2 <f> mn - 1) 



sin kRn 

(kR HI II 



cos kR m , 

(kRrnn)*' 



+ - (1 - cos 2 mn ) — — , (8) 

L Kii mn 



III. LARGE DETUNING LIMIT 

For red detuning (Aa < 0) the atoms are attracted to 
high-intensity regions of the field. In the large detuning 
limit, i.e., where the magnitude of A a exceeds the atomic 
linewidth 7 by far, | A^| 3> 7, the laser field creates a con- 
servative potential for the atoms. The recoil noise is so 
strongly suppressed that heating plays no role on the rel- 
evant timescale of motion. It is enough to consider only 
the conservative part of the dynamics described by the 
Hamiltonian. The dissipative processes will be taken into 
account later by the introduction of a phenomenological 
temperature. 

The internal electronic dynamics of the atoms consists 
of fast oscillations on a short timescale. This can be adi- 
abatically eliminated to derive its effect on the external 
motion. The adiabatic atomic polarizations, according 
to Eq. ([3]), obey the implicit equation 

Q N 

On = -7—f(r n ) + -7— ^ {Oimn ~ iPmn) <7 m . (12) 

Aa Aa 

To leading order in |7/A^|, the polarization is just the 
first term, which gives the zeroth order of the Hamilto- 



where <j) mn = Z(R mn ,x) is the angle between the dis- 
tance vector of the two atoms and the axis of the po- 
larization, x. As this function is singular at the origin, 
kRj n n * 0, later we will use rather its Fourier transform 



(3{k) = lim 



3 (k 2 z +k 2 )(k 2 -k 2 +r, 2 ) 
v^fc? (k 2 - k 2 ) 2 + 27f(fc 2 + k\) 



(9) 



B. Circular polarization 

If the polarization is circular, e — — -^(x + iy), the 
atomic quantization axis is the field propagation direc- 



H, 



(0) 



ctr 



N 

E 



-*2 

P£_ 

2m 



Ml 2 



hjn 2 
a 2 

A 



N 

E 

n , m 



/3(r n - r m )f(f n )f(fr, 



(13) 



This Hamiltonian H^g describes the dynamics of dipo- 
lar particles, where the polarization is induced by the 
external driving field. The neglected higher order terms 
describe the effect of one polarized particle on the polar- 
ization of another, i.e., the local field effects. The same 
effect is at the heart of the Lorenz-Lorentz refractive in- 
dex of a dielectric medium. We will later discuss in which 
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parameter regime the local field effects become significant 
and present a systematic derivation of an effective two- 
body Hamiltonian. 



IV. MEAN-FIELD APPROXIMATION 

In the following we will resort to a single-atom analysis 
based on a mean-field potential. The conservative dipole 
trap potential is 



The mean-field dipolc-dipolc potential is 



(14) 



r(0). 



o; 2 



-2^7^/(7=) / d s r 2 p(r 2 )l3(r-r 2 )f(r 2 ), 

(15) 

where the continuous position distribution p(f) was in- 
troduced. This distribution is normalized to the number 
of atoms in the trap N. The convolution integral can be 
evaluated as a product in Fourier space [23, [30] • 

The final step in setting up the model consists of 
assuming that the cloud of atoms is described by a 
canonical ensemble at an equilibrium temperature T. 
By this approximation we avoid describing the heating 
and cooling processes which lead to the steady-state. 
Both are slow processes as the recoil noise is strongly 
suppressed. The canonical distribution provides a self- 
consistent equation for the spatial distribution of the 
atoms: 



p(r) = exp 



V(f,p(r)) 



(16) 



where the partition function Z ensures that the integral 
of p(x) gives the number N of atoms. This self-consistent 
equation for the atom distribution is the basic equation 
of our model, which can be solved only numerically. The 
same method has been used to study phase transitions of 
atom gases when the atom-atom interaction is mediated 
by a cavity field [Mil]. 

The MF model can be expressed in terms of dimen- 
sionless parameters, which amounts to the identification 
of the relevant quantities describing the equilibrium of 
a trapped cloud of atoms. The dipole trap depth (see 
Eq. (fT4| ) is set by the intensity and the detuning, how- 
ever, in the self-consistent equation (| 1 6[) it is compared to 
the temperature. None of the above quantities appears 
separately, thus it is appropriate to introduce the scaled 
intensity 



X : 



nn 2 



\A A \k B T 



(17) 



The depth of the MF dipole-dipole potential relative to 
the trap depth is determined by the product of the small 
parameter j/Aa and the atom number N through the 



atomic density p(f) in Eq. (|I5p . The appropriate param- 
eter is then the scaled atom number 



TV = 



N-y 

Ja7\ 



(18) 



With these two dimensionless parameters, X and A/", the 
effect of all the relevant physical quantities can be de- 
scribed. 



STABILITY ANALYSIS 



The solution of Eq. p^|) can be determined numerically 
by iterating the atomic distribution as follows. Initially 
we take the canonical distribution of a noninteracting gas 
in the dipole trap. The dipole-dipole interaction term 
given by the convolution integral in Eq. (|f 5p is calcu- 
lated then from this distribution. It is suitable to go into 
Fourier-space where the convolution is a simple product. 
The Fourier transform of the term f3(r) is known, cf. Eqs. 
([9]) and (fTTj) . while the term p(r)f(r) is transformed by 
numerical Fast Fourier Transform (FFT). Due to the fi- 
nite support of the distribution p(k) and the mode func- 
tion f(k), the singularity of f3(r), which appears as a 
non-decaying Fourier transform in momentum space, is 
automatically regularized. On transforming the product 
back to real space by FFT we get the dipole-dipole term. 
Adding it to the dipole trap term of Eq. (|14[) yields the 
total mean field potential which furnishes a new atomic 
distribution via the canonical form in Eq. (| 16[) . The re- 
sulting p(x) can be used as the starting distribution in 
the next step of the iteration. Continuing the steps of it- 
eration until convergence, one obtains the self-consistent 
solution of Eq. (fTBj) . 

The iteration does not necessarily converge: instability 
can be induced by large enough atom number or intensity. 
The iteration method suggests that the instability occurs 
as a collapse of the atomic cloud due to self-contraction in 
the center of the trap. However, in the lack of a precise 
modeling of the collisional processes, the collapse itself 
cannot be accounted for by our approach. We must limit 
ourselves to determining the range of convergence. 

The mean field dipole-dipole energy Eq. (p~5|) depends 
on the shape of the atomic cloud. It is the pancake- 
shaped trap where the dipole-dipole contribution to the 
MF potential is negative, deepening the trap depth in the 
center [26j , so that the collapse of the cloud can be ex- 
pected. It vanishes in the center for a spherical atom dis- 
tribution (the refinement of this statement can be found 
in [371]). For a cigar-shaped cloud the MF dipole-dipole 
energy would be positive, repelling atoms from the cen- 
ter, and thus the instability we discuss in the following 
could not occur. 

In Fig. [T] the stability range of the iterations is pre- 
sented for circular polarization in the x-y plane, and for 
beam waists w = 1.33, 2.66, 5.33, and 10.66 in units of A. 
On the two-dimensional plot for the scaled atom number 
and scaled intensity variables ( "phase diagram" ) , border 
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FIG. 1: The boundary of the stability range of the atomic 
cloud is shown in logarithmic phase diagram of the scaled 
intensity I and the scaled atom number N '. The field polar- 
ization is circular in the x — y plane, and the waist is w/X= 
1.33 (circle), 2.66 (triangle), 5.33 (diamond), and 10.66 (pen- 
tagon). Straight lines represent a fit on the numerical data. 
As shown with empty circles, below T = 10 the phase bound- 
ary deviates from the corresponding fitted line (w/X = 1.33), 
which indicates the effect of untrapped atoms. 



FIG. 2: Left scale: Maximum density of the cloud, p(0) in 
units of 1/A 3 x Aa/j, plotted against the scaled intensity T 
for w/X— 1.33 (circle), 2.66 (triangle), 5.33 (diamond). The 
straight line is a fit on the w/X = 5.33 data. The scaled 
atom number M is set such that the system remains slightly 
below the critical point. Right scale: stars with a constant fit 
represent the MF dipole-dipole potential Vdd/ ksT calculated 
close to the boundary of the range of stability. The beam 
waist is w/X = 5.33. 



points of the stability region are shown and the conver- 
gent iterative solution of Eq. (|16p exists in the region 
below the points. The border points can be well fitted by 
a power law dependence of the critical scaled intensity 
on the scaled atom number, X oc J\f~ c . The fit is repre- 
sented by lines, the exponent is c = 0.40(±0.01). Shown 
only for the w/X = 1.33 data with empty circles, the 
boundary bends away from the fitted straight line at the 
right-most end. This happens below a certain scaled in- 
tensity (X < 10), when a significant portion of untrapped 
atoms appear. It is not shown on the figure, however, we 
obtained the same results for linear polarization along x. 
We also checked that the addition of an arbitrary con- 
stant to the Fourier transform {3(h) does not appreciably 
shift the phase boundary. This justifies the neglect of 
any type of Dirac-i5 (r) potential in the Hamiltonian, e.g., 
the contact potential [35| or s-wave scattering. 

The dipole-dipole coupling enhances the trap depth in 
the center and increases the atom density there. The re- 
sulting self-contraction of the cloud is counteracted by 
the random motion of the atoms. We expect that in- 
stability occurs when the energy shift due to the dipole- 
dipole coupling exceeds the thermal energy of the atoms. 
The numerical approach has allowed us to confirm this 
expectation. In Fig.[2]the ratio of the dipole-dipole inter- 
action potential and the thermal energy is plotted (with 
stars) on a logarithmic scale at the edge of the stable 
region (where still stable solution exists), and one finds 
Vdd ~ ksT closely constant. 

Accepting the instability condition Vdd ~ k B T, the 
critical exponent c = 0.4 can be obtained by simple ar- 
guments. It follows from Eq. QU that V^\o)/k B T is 
proportional to XM times the convolution integral. Us- 
ing the dominant term of the potential for the distribu- 
tion, i.e., the trap potential in harmonic approximation, 



the resulting Gaussian distribution has a normalization 
factor proportional to X 3 / 2 . The remaining part is the 
convolution involving the function (3 which is singular 
at the origin. Thus the main contribution must come 
from this small domain, the cloud size is irrelevant and 
the integral must be determined, at least to leading or- 
der, by the aspect ratio of the trapped cloud. Altogether 
V dd (°) l k B T « I 5/2 AA from which the scaling X oc J\f~ 0A 
follows. 

For consistency, let us check the validity of the effec- 
tive Hamiltonian given in Eq. (fT3"|) which keeps only the 
leading order term of the dipole-dipole interaction. Fig- 
ure [3] depicts the ratio of the MF dipole-dipole potential 
to the trap potential at the center of the cloud, which 
is the same ratio as that of the higher order terms of 
the polarization to the leading order one [cf. Eq. (fT2|) ]. 
This ratio grows as a function of the scaled intensity un- 




FIG. 3: The ratio of the MF dipole potential and the total 
potential at the origin as a function of the scaled intensity I 
for w/X = 5.33, and at fixed values of AF = 1, 2, 4 X 10~ 5 . 

til reaching the critical point which is clearly manifested 
in this semi- logarithmic plot. The main thing to observe 
is that the criticality is reached at fairly low ratio of the 
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dipole-dipole potential to the trapping one, at a value 
well below 0.1. Therefore, in the considered parameter 
regime, the leading order description is enough to find 
the boundary of the stability domain. For other parame- 
ter regimes, where the ratio would be larger than 0.1 but 
still below 1, we present, in the Appendix, a Hamiltonian 
which accounts for the two-body interaction to all order. 

VI. POSSIBILITY OF BOSE-EINSTEIN 
CONDENSATION 

The peak density is a key quantity with respect to 
reaching quantum degeneracy, which happens roughly at 
the density p(0)A| eBroglie « C(3/2) = 2.612 [H, or, ex- 
pressed in another way, 

p(0)A 3 «4l(^-) 3/2 , (19) 

where T lec = (h 2 k 2 /m) /fcs is the recoil temperature. 
This condition depends, separately, on the temperature 
which has been, so far, embedded in the scaled intensity. 
Therefore, depending on the working point, i.e. the actual 
value of the temperature and the trap depth, the order 
of reaching quantum degeneracy or reaching the dipole- 
dipole instability can be chosen. 

Obviously, both increasing the trap depth (at fixed 
temperature) and increasing the absolute atom number 
(at fixed detuning) lead to an increase of the density at 
the trap center. However, the critical behavior can be 
induced at different peak densities depending on how it 
is obtained: whether by increasing T or M . We present 
in Fig. [2] the calculated peak densities along the phase 
boundary, i.e., for clouds which are infinitesimally close 
to the boundary of the stable phase. The numerical val- 
ues obtained for various waists overlap, and can be well 
fitted by a power law function 1/X. Thus the dipole- 
dipole instability appears at densities 

P(0)A3 « . (20 ) 

7 ^ 

We recall that the dimensionless detuning A ,4/7 de- 
scribes the ratio of the dipole trap potential depth and 
the spontaneous photon scattering rate. Thus if the 
right-hand-side of the inequality ([2D]) is smaller than that 
of (fl"9)) . the dipole-dipole instability prevents that the 
limit of quantum degeneracy could be reached. For ex- 
ample, with atoms at about 100 /iK, 3 orders of mag- 
nitude above the (typical) recoil temperature, and for 
T = 10 (trap depth about mK), one needs a detun- 
ing \Aa\ > 10 5 7 so that the Bose-Einstein condensation 
point falls into the stability range of the dipole-dipole 
interaction. On increasing the scaled intensity T, the 
threshold density p(0) for the dipole-dipole instability 
decreases (see Fig. [2] and Eq. ([20])). which amounts to 
a more demanding condition on the atom number or on 
the temperature needed to achieve condensation. 



VII. CONCLUSION 

We set up a simple but versatile model to study the 
stability of far-off-resonance trapping against the dipole- 
dipole interaction. The model can be adapted to var- 
ious geometries, atomic species, polarizations, etc. We 
found that the radiative atom-atom interaction in the 
laser polarized gas can indeed produce an instability of 
a pancake-shaped atomic cloud. When the dipole-dipole 
energy shift reaches the temperature, the unbalanced at- 
traction in the trap center yields a collapse of the cloud. 
The condition for the instability in terms of the atom 
density and the temperature is obtained numerically in 
the form of phase diagrams. 

The density necessary for large enough dipole-dipole 
shift is quite large, however, the temperature is not 
needed to be very low so that the instability be observ- 
able. Provided the optical dipole trap is deep enough, 
the stability can be lost at moderate phase-space den- 
sities, well above the condensation threshold. On the 
other hand, it is also possible to find a working point, 
very large detuning, low temperature and not too large 
trap depth, where the dipole-dipole coupling itself does 
not prevent the gas from forming a quantum degenerate 
state. The dipole-dipole interaction strongly dominates 
the influence of quantum statistics in near-resonant laser 
fields This is not the case for very large detunings, 
hence we expect that the effect of cloud instability offers 
parameter regimes where quantum statistical properties 
can be probed. This is the subject of our forthcoming 
research. 
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VIII. APPENDIX 

In the Appendix we present an effective two-body 
Hamiltonian H c g that takes into account the possibility 
of close pairs of atoms in a non-perturbative manner. 

Eq. (fl"2)) can be solved iteratively to obtain 

A N 

= /(*n) + e «m/(r ro ) 

^ m— 1 
/ \2N N 

\ A ' m=l 1 = 1 
m^n l=£m 

, x 3 N N N 
\ A ' m=l 1=1 k=l 



7 



where r 



iflmn has been introduced for the effective two-body Hamiltonian: 



sake of compactness. The series includes terms of high 
order in the small parameter 7/|A^| <C 1, however, not 
all of them can be neglected. Such a high-order term 
can contain extremely large factors provided all the Qki 
in it have large values, i.e., the corresponding atoms all 
are close to each other (within a small fraction of the 
wavelength). This is because in the limit R m n - > 0, 
the dipole-dipolc coupling function gives f3(R mn ) — * oo. 
However, even in this limit a is bounded, i.e. a(R mn ) — > 
1, thus the real part of can safely be neglected in the 
large detuning limit. 

A systematic low-density approximation of Eq. (|21j) 
can be made by considering only two-body terms. This 
amounts to the assumption that any pair of atoms can be 
very closely spaced, however, in this case there is no third 
atom in their immediate vicinity. Thus if for any n =/= m, 
flnm ^> 1, then for every I = 1, . . . ,N such that I ^ n, 
I 7^ m, we have flnmPmi < 1- In this approximation, 
we reorder the sum, and obtain two geometric series of 
-l3 mn -f/A A , leading to 



AT 



1 + Et 



St 

m 



(J3 mnl /A A f 

- {i3 mnl /A A y 

((3 mnl /A A ) 



f(rn 



((3 mnl /A A y f{rn 



(22) 



Here we have assumed that for every atom pair n =/= m, 
/ A A \ < 1. This solution, within the MarkofT and 
the adiabatic approximations, takes into account in a 
non-perturbative manner the pairs of atoms separated 
by small distances. 

Inserting back the stationary adiabatic atomic polar- 
ization a n of Eq. ([22]) into the original form of the Hamil- 
tonian of the system Eq. ([3]) and eliminating double sums 
on the same ground as above, one can end up with the 



tt _ \^ Pn , nAl 



n,m 

E /m/nAnn { (l + P mn ) ' + 3(/? m „) 2 ) 



(23) 



where /„ = f(f n ) is the mode funcion of the field at the 
position of the nth atom, and /3 mn , (3 mn are defined as 
below: 



fa 
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A, 



A A i 



ft 



£2 h'mn 



(24) 
(25) 



Note that the denominator of the expressions Eq. (|24|) 
and Eq. (|25| should be strictly greater than 0, due to the 
relation (3 mn j/\A A \ < 1 that is the condition of summing 
up the geometric series resulting from Eq. (|21[) . 

Apart from the zeroth order effective Hamiltonian 
•^iff > Eq. (fl3|) , further terms appear in Eq. ([23]) in higher 
orders of (3 m n) 'I ' A A . The new terms have two appar- 
ent meaning: the second line of Eq. (|23|) renormalizcs 
the dipole trap potential (first line), while the appear- 
ing j3 mn , f3 mn in the third line increase the strength of 
the dipole-dipolc interaction. This increment is negligible 
and the zeroth order effective Hamiltonian, Eq. (flU|) , can 
be used if [3 mn "// A A is really a small parameter, however 
when it gets close to one, the non-perturbative form of 
the effective Hamiltonian, Eq. (|2"3")) , should be used. 
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